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Heavy quarks and the collective properties of hot QCD 
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After reviewing the evidence that QCD matter at ultrarelativistic energies behaves as a very good fluid, we 
describe the connection of QCD fluidity to heavy quark observables. We review the way in which heavy quark 
spectra can place tighter limits on the viscosity of QCD matter. Finally, we show that correlations between flow 
observables and the event-by-event charm quark abundance ("flavoring") can shed light on the system's equation 
of state [1]. 



The purpose of high energy heavy ion colli- 
sions is to study the thermodynamic properties 
of strongly interacting matter at high tempera- 
tures and densities. The most ambitious part of 
this program is to create a bubble of deconfined 
"quark-gluon plasma" , a gas of "free" quarks and 
gluons mimicking the properties of the universe 
shortly after the big bang 

This program, however, has a potential ob- 
stacle: Thermodynamics generally applies to 
"large" ," static" systems, where "large" is com- 
pared to typical sizes of microscopic degrees of 
freedom. While the system created in a heavy 
ion is large compared to the quark or hadron size, 
the system is certainly not static. In fact, it could 
be be a far-from-equilibrium "mess of quarks and 
gluons" , where each parton scatters a few times, 
but whose dynamics has no connection to ther- 
modynamics. In this case, concepts like phase 
transition have little value. 

The best that we can hope for is to replace 
"static" by "slowly evolving", where "slowly" is 
defined in respect to the timescales of the mi- 
croscopic processes within the system. If the 
system is slowly evolving (alternatively, if mi- 
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croscopic dynamics equilibrates very fast, as it 
naively should in a strongly coupled hot sys- 
tem), we can explore the thermodynamics of the 
system via hydrodynamics:The systems collective 
motion will be governed by the equation of state 
and transport coefficients, calculated from equi- 
librium thermodynamics. 

RHIC experimental data has given us reason 
for optimism in this respect: One of the most 
widely cited findings in ultrarelativistic heavy 
ion collisions concerns the discovery of a "per- 
fect fluid" in collisions at the Relativistic Heavy 
Ion Collider (RHIC) [SI]. The evidence for these 
claims comes from the successful modeling of the 
anisotropic expansion of the matter in the early 
stage of the reaction by means of ideal hydrody- 
namics |5l6l7j 

It appears, therefore, that the system is, to 
a very good approximation, locally thermalized, 
and we can use flow observables as a probe of the 
equation of state, and eventually of phase transi- 
tions. However, efforts in this direction have met 
remarkably little success: The probes for which 
hydrodynamics works very well are insensitive to 
the equation of state [617] . Probes (such as HBT 
radii and average transverse momentum) which 
are more sensitive show no structure indicative of 
a phase transition |5I6I7| , presumably because the 
systems initial conditions are not "tidy" enough 
to highlight discontinuities in the equation of 
state 
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Other probes more sensitive to the details of 
the equation of state are therefore needed. One 
obvious candidate is heavy quark observables. A 
heavy probe in a thermahzed medium is, as is 
well-known, described by Brownian motion. This 
makes it a very sensitive probe of thermalization, 
since the scale of thermalization of a heavy parti- 
cle Theavy IS paramctrically larger wrt to its light 
component Tught ^ v/iTs) (where ij/s is the di- 
mensionless ratio of viscosity to entropy density 
and T is the temperature): If the mass of the 
heavy particle is M > T, Theavy / night ^ M/T 

m- 

RHIC data to date is consistent with "heavy 
particles flowing as much as light particles" [11] 
[12], although large systematic errors persist (for 
example, we can not distinguish a charm from 
a bottom quark) If this is really true, than rj/s 
is truly <C 1 and heavy quarks are, for all in- 
tents and purposes, thermahzed with respect to 
the rest of the system. They (in particular, the 
more abundant charm quarks) can then be used 
as chemical probes for the equation of state. 

In general, charm in heavy ion collisions is not 
expected to be chemically equilibrated. The bulk 
of charm content should be produced by "hard" 
processes in the initial state at a concentration 
far above their equilibrium expectation [13] . The 
abundance of cc pairs produced in heavy ion 
collisions at the LHC is expected to be reason- 
ably high [13] (~ 10^~^, of course paramctrically 
smaller than the total ^ 10* multiplicity). In the 
dilute limit, the total charm abundance is for all 
intents and purposes a conserved number. Hence, 
a good observable is the dimensionless quantity 
p = p/s ^ pT~^ <^ 1. The events entropy, also 
nearly conserved for a good fluid, can be related 
to the multiplicity rapidity density dN/dy [2]. 

In the dilute (corrections of O (p^)) infinitely 
heavy quark (corrections of O {T/M)) limit, the 
contribution of an abundance of charm quarks 
to the free energy can be computed by adding 
a Polyakov loop density [T4l to the free energy 
density 

J^{T) = MT) + PSo{T)Fq{T) (1) 

where Fq{T) = -T\n£{T) and e{T) is the renor- 
malized Polyakov loop. The Polyakov loop is ob- 
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Figure 1. (color online) Top panel: The de- 
pendence temperature, for p = (dashed lines) 
and p = 0.1 (solid lines). The width of the bands 
denotes lattice uncertainties [15] in the speed of 
sound of a 2-1-1 QGP, while the thick line de- 
notes the conformal, non-interacting value where 
cl = 1/3. Bottom panel: p vs -/s for Pb-Pb col- 
lisions computed using Eq. dS]). The dashed lines 
denote the uncertainties in the perturbative QCD 
calculation 



tainable from lattice calculations, and the quan- 
tities underscored with denote the values be- 
fore the charm flavor was included. The speed of 
sound can then be computed by textbook thermo- 
dynamic methods, cl — j^f^ , and s — —dF/dT. 
Fig. [T]top panel shows our estimate for the speed 
of sound derived via Eq. ([1} using the expectation 
value of the Polyakov loop extracted from the lat- 
tice (2-1-1 QGP with almost physical quark masses 
[15]). One can see that the main effect comes 
from the region near the phase transition (where 
there is a minimum in the speed of sound) but 
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well before the Polyakov loop expectation value 
reaches its asymptotic high-T limit, leading to a 
negative shift of the speed of sound from its value 
in a 2+1 QGP. This can be physically readily un- 
derstood: Correlations between the medium and 
slowly moving heavy quarks lower the system's 
response to pressure. 

Our estimate stops close to Tc, as —T\n£{T) — > 
oo in the confining phase (where the Polyakov 
loop expectation value vanishes). Mathemati- 
cally, one can trust our approach as long as the 
heavy quark is much heavier than any other scale 
in the system, ie -Tln^(T) < M,. At some 
point in the approach to confinement, however, 
this approximation breaks down. To estimate the 
contribution of flavoring in the confined phase we 
assume flavorful confined QCD is described by 
the hadron resonance gas model. In this case, fla- 
voring can be approximated by an admixture of 
heavy mesons in a gas of pions. The latter has a 
speed of sound of c'/"'*' ~ l/V^ (ultra-relativistic 
ideal gas), while the former will have a speed of 
sound of c'^'^'^^y ~ ^/5T/{3Mmeson). The speed 
of sound of the mixture will therefore go as 

^'-l + O (P^f^) (2) 

parametrically smaller than the contribution in 
the deconfined phase, which is just O (p). Thus, 
the flavoring effect on the speed of sound is spe- 
cific to the deconfined phase. The effect of fla- 
voring in a weakly coupled QGP is similarly de- 
scribed by Eq. [2] (with M„eson Mcharm), so a 
correlation of Cs with flavoring would indicate a 
deconfined but strongly coupled system. 

These effects produce observable consequences. 
Using pQCD estimates TJ" and logarithmic scal- 
ing for multiplicity rapidity density, we estimate 
p at the LHC to be (Fig. [1] bottom panel) 

, _ 1 dN charm /dy 
6 dNcharged/dy 

1 NcoUisions '^pp-i-cc 

qTT 7 ? ^ ^-^^ 

<J -^^participants Yn 

This is the average expectation. dNcharm/dy and 
dN/dy will however vary event by event, in an 



approximately Poissonian manner. Provided that 
charm can be reasonably reconstructed and there 
is a large enough event sample, p is an experi- 
mental observable capable of serving as a binning 
class for events (see Fig. [2]). 

As is well known, there is a connection be- 
tween the speed of sound and the limiting aver- 
age velocity of a hydrodynamic expansion with 
shock-like initial conditions, {itVt) freezeout ~ 
f{Npart) (cs)t- where "freezeout" implies averag- 
ing over the freeze-out hypersurface [2] while the 
subscript r means the average is done over the hy- 
drodynamic evolution. For a shallow shock this 
result is exact [TB]. While knowledge of the ini- 
tial geometry is needed to establish the form of 
f{Npart), model calculations [6| indicate that the 
dependence is not washed away even in steeper 
shocks and more complicated initial geometries. 

The final transverse flow is in return connected 
to the average transverse momentum (px) — 
T + ni (jtVt)- Hence, the decrease of the speed 
of sound close to Tc (Fig. [1]) could lower (px) 
for more flavored events with respect to flavor- 
less ones (Fig. [2] right panel). Note that this 
effect is opposite to the naturally expected pos- 
itive correlation due to the correlation between 
NcoUisions, N part and dN/dy. The coefficient as- 
sociated with this heavy flavoring effect would 
be straightforwardly related to non-perturbative 
QCD via Fig. [T] This effect might be easier to 
measure in smaller systems due to the greater 
event-by-event variation in p and less background. 
The main requirement of such an analysis is the 
ability to experimentally gauge both the charm 
quark abundance and (pt) event-by- event 

A possible "trivial" effect which would give 
correlations in the same sense as the effect pro- 
posed here is energy conservation (roughly, charm 
quarks need a lot of energy to be created, and that 
lowers (pt) ) ■ The correlation due to energy con- 
servation should, however, be suppressed by fac- 
torization and boost-invariance: Charm quarks 
are created of partons with larger Bjorken x jl7j 
than the mid-rapidity soft particles accounting for 
most of dN/dy, so the energy they take up comes 
from regions where rapidity deviates from zero. 
We do not expect, therefore, that energy conser- 
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Figure 2. (color online) Left panel: The event-by-event distribution of events with respect to charm 
content and Npart, together with the cut required to analyze the response of the system to charm. Right 
panel: The dependence of (pt) with charm number. 



vation will lower (pr) at mid-rapidity. Hence, the 
observation of a significant charm- flow correlation 
should be due to the collective properties of the 
system, not energy conservation. 

In conclusion, we have argued that heavy quark 
observables seem to confirm that the matter pro- 
duced in RHIC collisions is a very good fiuid. 
This means that heavy quarks are, too a good ap- 
proximation, in thermal (but not chemical) equi- 
librium with the rest of the medium. This makes 
it possible to study the medium response to the 
presence of heavy quark impurities. This re- 
sponse is calculable, to a good approximation, 
from lattice QCD data, and can produce non- 
trivial correlations between event charm abun- 
dance and flow properties. 
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